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Abstract 



We study a motion of quantum particles, whose properties depend on one 
coordinate so that they can move freely in the perpendicular direction. A 
rotationally-symmetric Hamiltonian is derived and applied to study a gen- 
eral interface formed between two semiconductors. We predict a new type 
of electron states, localized at the interface. They appear whenever the two 
bulk dispersions intersect. These shallow states lie near the point of inter- 
section and are restricted to a finite range of perpendicular momentum. The 
scattering of carriers by the interface is discussed. 
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I. INTRODUCTION 



Highly developed methods of crystal growth allow building of mesoscopic systems with 
spatially varying parameters. The motion of carriers in such a system is described within 
the effective mass theory (EMT) by a quantum equation. Depending on the semiconductors 
chosen, this can be either Schrodinger-like for wide-gap materials, or Dirac-like or more 
complicated for narrow-gap semiconductorsQ. Thus, there is a way to realize in practice 
the quantum mechanics with parameters depending on coordinates. All information about 
the periodic quickly varying potential of the lattice is contained in the effective mass or the 
"velocity of light" - the interband matrix element. 

The foundation of the phenomenological EMT is based on the envelope function approx- 
imation (EFA). The total wave function is represented as the product of the Bloch function 
and a slowly varying envelope function; it is the latter which enters the quantum equation. 
The commonly used heuristic efa! exploits the fact that the interband matrix element 
does not change much in many important for practice III-V semiconductors. This leads to 
a reasonable assumption that the Bloch function can be chosen to be identical through the 
whole system. The resulting EMT is equivalent to the conventional quantum mechanics, e.g. 
the envelope function should be continuous at any heterojunction. However, the interband 
matrix element does vary for different semiconductorsi. Recently, a new exact version of 
the EFA was suge estecfl, where only the total wave function has to be continuous and the 
envelope function does not need necessarily to fulfill this requirement. 

In Ref.i, first evidence indicated the existence of new interface states, impossible in the 
conventional quantum mechanic,. The further investigation! revealed a novel effect for the 
coordinate-dependent two-band Dirac Hamiltonian: A localization can occur at a step-like 
potential (junction), provided that the free motion along the junction is accounted for. The 
envelope wave function was shown to be discontinuous for these states. In the present paper, 
we study an arbitrary junction of two wide-gap semiconductors. Also here, the free motion 
along the junction leads to binding in its perpendicular direction. 

Although graded crystals with very different dependences of the effective mass m(z) on 
the coordinate z can be grown, the most attention has been drawn to abrupt heterostruc- 
tures. This is probably because one does not need to use a Hamiltonian with coordinate- 
dependent kinetic term (see below). Instead, a problem of choosing the matching conditions 
(MC), connecting the wave functions on either side of a junction, is met. Usually they are 
derived from the requirements of the continuity of the probability and its current at the 
junction. These requirements, although being correct, are not constructive. Indeed, they 
are meaningless for matching of exponentially growing or decaying solutions as their current 
vanishes. On the other hand, incident and scattered waves appear in the propagating case. 
Two continuity equations are not enough to fix the MC and two constants of the scattering 
problem. 

An alternative approach of deriving the MC directly from the Hamiltonian is suggested 
in Ref.i: The abrupt junction can be treated as a limiting case of a smooth one, when its 
"smoothness" length tends to zero. For the graded case, the Hermitian Dirac operator is 
derived and from it the MC are extracted. Apart of appealing to intuition, such a method 
is shown to be equivalent to the conventional scheme. The Dirac Hamiltonian is conve- 
nient to handle with, because it contains first-order derivatives, so that it can be easily 
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symmetrized and it is originally rotationally-symmetric. In the present paper, we follow the 
above- described approach and develop an effective way to derive Hermitian, rotationally- 
invariant one-band Schrodinger Hamiltonian. The MC depend on the free motion, since it 
is incorporated in the Hamiltonian. 

The paper is organized as follows: In Sec. II, we derive the Hamiltonian and obtain from 
it the matching conditions. In Sec. Ill, the dispersion of the interface states localized at the 
hetero junction is found and the scattering of the carriers by the junction is discussed. Sec. 
IV summarizes the results. 

II. ROTATIONALLY-SYMMETRIC SCHRODINGER HAMILTONIAN 

We study a junction of two semiconductors with the effective masses rrii and the bandgaps 
2Aj, dependent on one coordinate z, i = 1 for z < and i = 2 for z > (see Fig. 1). Since 
all parameters depend on z only, the free motion in the perpendicular direction x has a 
good quantum number k±_ so that the wave function is ip ~ exp(ik±x). (Hereinafter we put 
h = 1.) Let us start with the case of vanishing free motion k± = 0; we consider k± ^ later. 
The problem one immediately meets while studying a graded crystal with some function 
m(z) is that the kinetic term T = d 2 /2m is not Hermitian. Its most general form, after 
symmetrization, was foundS to equal 



with constants a and (3 such that 2a + (3 = —1. Note the appearance of singular terms 
for abrupt junctions m" ~ S'(z)A m and (m') 2 ~ S 2 (z)(A m ) 2 , where S(z) is the Dirac delta- 
function and A m = m 2 — These singularities lead to a discontinuity of the wave function. 
Representing the wave function if)(z) = ip(z)f(z) as a product of the continuous part /(+0) = 
/(— 0) and the function (f(z), which jumps at the junction, A^ = (p(+0) — y?(— 0), one 
obtains ip' ~ fA^S(z) and ip" ~ fA v 5'(z) + 2A tp S(z)f. Substituting these expressions 
into Eq. (|2.1|) and selecting the singular contributions, proportional to 5'{z) and S 2 (z), one 
finds that all the singularities in Eq. ( |2.1| ) are compensated under the following conditions: 
A,p/(p + aA m /m = and a{3(A m )/m — A v /ip = 0. These give either (3 = — 1 or a = 0. 
Recalling that 2a + (3 = —1, one notices that both possibilities lead to the same effective 
Hamiltonian T(z) = d z (l/2m\dz. acting on the continuous function f(z). This form was 
already derived in many worksHll 2 ], where arguments different from ours were used. 

The next step is to incorporate the free motion along the surface. Up to our knowl- 
edge, the Hermitian rotationally-symmetric form of the Schrodinger Hamiltonian has not 
been written yet. The only form, being simplified for k± = to the one- dimensional case, 
recovering the Laplace operator for constant m, and invariant with respect to orthogonal 
rotations involving z and x is given by 



where A = ±1. The physical meaning of the value A will be elucidated later, from the 
comparison of the one-band and the two-band Hamiltonians. The appearance of A is due 



T = m a d z m^d z m a /2 = — d 2 z - — d z + af3 ( — ) + a— , 

2m m \m J m 



(2.1) 



T=(d 2 -i\d x )^-(d z + iXd x ) , 



(2.2) 
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to the fact that the axial symmetry possessed by the hetero junction is broken as soon as 
k± ^ 0. Two signs of A formally account for two possible ways of writing the kinetic form 
symmetrically. Acting by the operator T on the function / = g(z) exp(ikj_x), we derive the 
eigenvalue equation for the energy e and half the energy gap, A(z), playing the role of the 
potential at k± = 0: 

( J_# - ^-d z + X^k ± -*L-A + e)g = . (2.3) 
\2m z 2m 2 2m 2 2m / V 1 

The presence of the third term in Eq. (|2.3|) , originating from the mixed derivative d z (l/m)d x 
in Eq. ( |2.2| ), will be crucial for the following discussion. Let us search for solutions of 
Eq. ( |2.3|) in the form g = exp[/ dz k(z)}. Far from the junction, z — > ±oo, the sec- 
ond and the third terms of Eq. fl2.3| ) vanish. The remaining terms determine the values 
Ki = ±^/2mj(Aj — e) + kj_, nothing else but the (imaginary) wave numbers of the two bulk 
semiconductors. The terms n' /2m — urn' /2m 2 + Xk±m' /2m 2 prevail at the junction. To avoid 
unphysical discontinuities in the abrupt case, we require [(k — Xk±)/m\ = 0, recovering thus 
the matching conditions 

[(k — Xk±)/m] l = [(— k — Xk±)/m\ 2 . (2.4) 

The sign of K\ and k 2 is chosen in the way to ensure a localized wave function. For the case 
m = const, these matching conditions coincide with those of BastardE! and do not depend 
on k±. 



III. INTERFACE STATES 

If the energy e is measured from (Ai + A 2 )/2, one rewrites k 2 = k\ — 2mi(e ± A), where 
the sign "+" corresponds to i = 1, "-" to i = 2 and A = (Ai - A 2 )/2. Equation Q, 
together with the bulk values of yields 

e(k ± ) = -P( mi + m 2 ) + X2k ± VP , (3.1) 

where P = |A/(m 2 — m\)\. To obtain Eq. ( |3.1| ) one squares Eq. (|2.4|), therefore it should 
be found which branch of Eq. (|3.1|) in which region represents the solution of Eq. (|2.4|) 
and corresponds to localized interface states. A simple analysis reveals that the states 
are localized in the region &;™ m < k± < fc™ ax , between minimal and maximal values of 
k 2 ± = AmfP. At k i± , the curve e(k±) is tangential to the bulk dispersions = ^A + k 2 _/2m i . 
The necessary condition for the interface states to exist is (mi — m 2 )(Ai — A 2 ) > 0, i.e. the 
bulk dispersions must intersect. The dispersion e(k±), Eq. (|3.1|) , and the bulk dispersions e« 
are shown in Fig. 2 for the following parameters: Ai = 1.2 eV, A 2 = 1.4 eV, mi = 0.01 m , 
m 2 = 0.02 mo (mo is the free mass of the electron). 

It follows from Eq. (^.4[) that Xk^(m 2 — mi) is positive. Recalling that the signs of A "+" 



and "— " originate from two possible ways of writing the Hamiltonian ( [2.2|) , we conclude 
that the opposite signs in Eq. ( |2.4j ) correspond to localized states for positive and negative 
k±. That is, if m 2 > mi, then k± > is localized for A = +1 in Eq. (|2.2j ) and k± < is 
localized for A = —1. The signs in Eq. ( |3.1| ) have to be chosen accordingly. 
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Since both values of A = ±1 are equivalent, the degeneracy of bulk dispersion of the 
Schrodinger Hamiltonian with respect to the sign of k± is preserved for the localized states: 
e (^_i_) — e (~ k±). This differs for the two-band model with the non-relativistic Dirac 
HamiltoniarJll, where the eigenvalues of energy e are classified (apart of momentum k±) 
by eigenvalues of helicity A = ±1. Here, the initial degeneracy of the bulk dispersion is 
lifted for the interface states: €\(k±) = €-\(—k±$L It is not difficult to check that the 
Klein-Gordon equation, obtained in the two-band model, contains the kinetic term analo- 
gous to Eq. ( |2.2|) . Thus, the value A in Eq. fl2.2|) corresponds to helicity in the two-band 
model. In discussing further similarities of the one- and two-band description, it is worth 
noticing following: Although qualitatively the dispersion of the interface states is similar 
in both models, corresponding to localized states between tangency points with the two 
bulk dispersions, the dispersion in the two-band model is nonlinear due to different bulk 
spectrum. 

Let us now briefly discuss the scattering of carriers by the interface. Like for the bound 
states, it is convenient to represent the wave function in the form g(z) = exp[±i / dz k(z)}. 
Then for z < it will be the sum of incident and reflected waves g = exp[i / dz Ki(z)] + 
B exp[—i J dz K\(z)] and for z > the transmitted wave is g — Aexp[i J dz K 2 (z)), where 
k\ = Ivriiie ± A) — k\. The condition of continuity of the wave function g yields 1 + B — A 
and the analog of Eq. (|2.4|) reads [(k± + in)/m + B(k± —iK)/m]\x = A[(k± + in)/m}\ 2 . The 
transmission coefficient D(e) = \A\ 2 k, 2 /k,i is given by 

D(e) = AK 1 K 2 m\/[k\(m 1 - m 2 ) 2 + {Kim 2 + K 2 mi) 2 } . (3.2) 

Far away from the point of intersection, k\ <C km\m 2 P or k\ 3> 4mim 2 P, function 
D(e) vanishes as square root as the energy tends to its threshold value (one of the bulk 
dispersions): D(e) ~ Kj. On the other hand, near the point of intersection, k\ m 4mim 2 P, 
the function D(e) vanishes linearly: D(e) ~ e — (mi + m 2 )P. This behaviour is typical for 
the conventional quantum mechanics: the coefficient of penetration D(e) is governed by the 
asymptotic behavior of the potential U(z) only. If the value U(+oo) — U(— oo) is zero than 
D{e) is linear at the penetration threshold. In the opposite case it vanishes as square rootLl. 



IV. CONCLUSIONS 



It is worth emphasizing that it is the free motion along the junction, which leads to the 
effects discussed. Although for k± = we have a step-like potential U(z) = A i; its role 
for k± is taken over by some more complicated function U ± (z,k±). Substituting into 
Eq. (|2.3|) g(z) — m 1 ^ 2 y{z) 1 we obtain a Schrodinger equation for the function y(z) with the 
potential 

C/ ± (^, kx) = 3m' 2 /8m 3 - m"/4m 2 - Xk ± m'/2m 2 + k 2 ± /2m + A . (4.1) 

The potential U form(z) = (mi+m 2 ) /2+(m 2 — mi) tan(z//)/2 with the smoothness param- 
eter I is presented in Fig. 3 as a function of z for several values of k± and I. The asymptotic 
behavior of the potential is preserved for any smooth m(z), C/(±oo) = A, + k 2 _/2mi. If 
the bulk dispersions intersect, then at k±, corresponding to the point of intersection, we 
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have no step in the potential any more: U(+oo) = U(—oo). As it is known from quantum 
mechanics, any potential with this property and containing a well, has at least one localized 
level. This level becomes shallow, if the potential changes rapidly, I — *> 0, (e.g. a 5-function) . 
For the abrupt potential, the part of the well is played by the matching conditions Eq. Q2.4|) . 
For k± 7^ and small Ki, i.e. near the intersection point of the bulk dispersions, Eq. (|2.4|) 
has a solution, describing a shallow level. The scattering by the junction is interpreted 
analogously. 

The states considered lie above the band edge, being embedded in the continuum spec- 
trum. However, they are real localized states, not resonant ones, since they lie in the energy 
gap of the whole system: Their energies are below the bulk dispersions, corresponding to the 
same value of k±. Naturally, many-particle effects, due to e.g. impurities and the boundary 
roughness could strongly affect these states. However, qualitatively it is clear that as long 
as these can be treated as elastic scatterer, the states preserve because in our analysis the 
absolute value of k± was essential, but not its direction. 

In summary, we have derived a rotationally-invariant effective one-band Hamiltonian, 
which mass depends on one coordinate. We applied it to derive the matching conditions 
for a single general junction of two semiconductors. They depend on the free motion along 
the junction, due to its mixing with the motion in the direction of growth. If the bulk 
dispersions of the two semiconductors intersect, then shallow localized states occur. The 
problem studied can be classified as an example of weak localization. The essential condition 
for the existence of the states, the intersection of the bulk dispersion curves, holds for a large 
variety of semiconductors!. Although more investigation is needed to understand in detail 
the many-particle effects, we believe that the states are not crucially sensitive to them. 
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FIGURE CAPTIONS 



Figure 1. Band-edge profile of an heterojunction. The energy gap is shaded. Because of 
symmetry, only the region of energy near the conduction bands is studied. 

Figure 2. Dispersion e(k±), Eq. (|3.1|) , of the system shown in Fig. 1 (solid curve; for 
parameters, see text). Dotted curves: bulk dispersions e\,2{k±)- Energy is measured from 
the middle of the bandgap. States are localized between tangency points ki± with the two 
bulk dispersions (heavy solid line). It is understood that similar curves occurs for k± < 
and near the valence band. 

Figure 3. Effective potential U + , Eq. Q4.1| ), of the smooth junction; for parameters: 1) 
k ± = 10 6 cm" 1 , / = 1.5 x 10 5 cm" 1 ; 2) k± = 5 x 10 6 cm" 1 , / = 1.5 x 10 5 cm" 1 ; 3) 
k ± = 5 x 10 6 cm" 1 , I = 10 5 cm" 1 ; 4) k± = 6.2 x 10 6 cm" 1 , U 1.5 x 10 5 cm" 1 . Note 
appearance of a barrier for A = —1, U~ (curve 5); parameters are the same as for curve 2. 
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Figure 1 
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Figure 3 
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